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We present a new approach to the induction detection of electron spin resonance (ESR) sig- 
nals exploiting the nonlinear properties of a superconducting resonator. Our experiments employ 
a yttrium barium copper oxide (YBCO) superconducting stripline microwave (MW) resonator in- 
tegrated with a microbridge. A strong nonlinear response of the resonator is thermally activated 
in the microbridge when exceeding a threshold in the injected MW power. The responsivity factor 
characterizing the ESR-induced change in the system's output signal is about 100 times larger when 
operating the resonator near the instability threshold, compared to the value obtained in the linear 
regime of operation. Preliminary experimental results, together with a theoretical model of this 
phenomenon are presented. Under appropriate conditions nonlinear induction detection of ESR can 
potentially improve upon the current capabilities of conventional linear induction detection ESR. 



Electron spin resonance (ESR) is a well-known method 
enabling direct measurements of the electron spin Hamil- 
tonian, with applications ranging from biology to mate- 
rials science and physics jlnSj]. However, a significant 
drawback of conventional ESR is its relatively low sen- 
sitivity compared with other spectroscopic and analytic 
techniques, such as fluorescence and mass spectrometry. 
For example, the world record in electron spin sensitivity 
stands today at ^ 10^ spins per 1 sec of acquisition (often 
denoted as spins/ V Hz) or slightly more than 10* spins 
in a reasonable ~ 1 h of acquisition which is still far 
from single electron spin sensitivity. This current sensi- 
tivity limitation also restricts the available imaging res- 
olution of heterogeneous samples. Thus, while the laws 
of physics do not set a limit to the spatial resolution of 
ESR (at least up to the atomic- length scale), in practice, 
as the image's voxel (volumetric pixel) size decreases, it 
contains less and less spins and thus quickly runs into 
the sensitivity limitation wall. For example, the systems 
in our laboratory, achieved recently a 440 nm resolution, 
limited mostly by spin sensitivity [5|. 

The above-mentioned numbers for sensitivity and res- 
olution refer to ESR systems that employ "induction de- 
tection", namely, they make use of Faraday's law for the 
detection of ESR signals by means of a pick-up coil or 
a microwave (MW) resonator. Induction detection is 
the basic principle behind all commercial state-of-the-art 
ESR systems; it enables the acquisition of high resolution 
spectroscopic data with complex pulse sequences; facili- 
tates the use of efficient imagining methodologies (mean- 
ing that signals are acquired and averaged in a parallel 
fashion from the entire sample); and features convenient 
sample handling. While our work is focused on induction 
detection ESR, other groups have looked into alternative 
detection methods in an attempt to increase sensitivity 
and resolution. These include, for example, magnetic 
resonance force microscopy scanning tunneling mi- 
croscopy ESR 01, spin-polarized STM electrically de- 
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tected magnetic resonance Q , and indirect spin detection 
via diamond nitrogen- vacancy (NV) centers While 
these new techniques may improve even more in the fu- 
ture, they have some inherent limitations, resulting in 
limited applicability. 

It is therefore evident that there is still a strong need 
to greatly improve the sensitivity of induction detection 
ESR up to the ultimate single spin sensitivity (in a rea- 
sonable acquisition time of ~ 1 h or less) , making it a gen- 
erally applicable method for noninvasive detection and 
imaging of small numbers of electron spins. Here we take 
a step in that direction and show that the sensitivity of 
induction detection may be enhanced if it is employed 
in conjunction with a unique new class of non linear su- 
perconducting resonators. In our new scheme, the sam- 
ple's resonance properties affect the non linear proper- 
ties of the resonator, thereby resulting in a complete new 
approach to the detection of ESR signals. Our initial 
experimental results demonstrate this new approach in 
practice and are accompanied by a theoretical analysis 
explaining our observations. The measured responsivity 
of our resonator is enhanced by a factor of up to 100 
when operating the system in the nonlinear regime. 

The experimental setup is schematically presented in 
Fig. [TJ It is composed of a stripline waveguide MW 
resonator [ll|, with a characteristic impedance of 50 fi. 
Here we employed its second mode at 6.1 GHz. The res- 
onator is weakly capacitively coupled to a feedline. In 
the resonator, a narrow section with a length of 11.5 jim 
and width of 0.3 pm is defined. This section is referred 
hereafter as the "microbridge". 

We began the fabrication process with a 0.5-mm-thick 
sapphire wafer coated with 150 nm of yttrium barium 
copper oxide (YBCO), and a gold cup layer of 200 nm. 
Gold contacts and alignment marks were patterned using 
standard photolithography and wet etching. The waveg- 
uide was patterned with electron beam lithography and 
wet etching. In the final step, the microbridge was pat- 
terned using a focused ion beam (FIB) system (T2| . 

The resonator was covered with another 0.5 mm thick 
sapphire wafer. It was then packed in a copper box whose 
inner top and bottom surfaces are covered with nio- 
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bium to provide superconducting ground planes for the 
stripline structure. Before installation, a l-mm-diameter 
hole was drilled on the top sapphire wafer, just above the 
microbridge. The hole was then filled with a common sta- 
ble free radical powder (DPPH from Sigma- Aldrich) . In 
addition a superconducting wire was placed between the 
two sapphire wafers, below the hole, and perpendicularly 
to the stripline. The wire allows applying low-frequency 
magnetic field modulation to the paramagnetic sample. 

The device was fully immersed in liquid helium inside 
the core of a superconducting coil. The ESR signal was 
measured in a continuous wave (CW) mode, similar to 
any conventional ESR system, with the exception that 
automatic frequency control was not used to track the 
resonator's resonance frequency. A monochromatic MW 
signal at a frequency close to the resonance frequency was 
injected. The reflected signal was measured while the 
external magnetic field was slowly scanned. In another 
setup, the resonator 6*11 reflection coefficient was mea- 
sured using a network analyzer. The resonator's loaded 
Q factor was found to be on the order of 5000. 

In previous work we have shown that the microbridge- 
integrated resonator system has a non linear response 
to a coherent tone injected close to the resonator's reso- 
nance frequency [l3| . When input power exceeds a cer- 
tain threshold, the power reflected off the resonator shows 
a self-excited amplitude modulation. We will refer to 
this threshold hereafter as the "modulation threshold". 
In Fig. [5] the spectral density of the reflected signal is 
demonstrated, above (panel b) and below (panel c) the 
modulation threshold. As will be shown below, the ESR 
responsivity can be enhanced when operating the res- 
onator close to the modulation threshold. 

We began by examining the linear response of the res- 
onator to the ESR signal. The resonator was excited with 
power well below the modulation threshold, while the 
static magnetic fleld was slowly scanned (O.OlT/min). 
We measured the reflection coefficient of the resonator as 




Figure 1. Experimental setup for ESR signal measurements. 
The device is installed in a cryostat where a strong static mag- 
netic field Bz is applied. A single coherent tone, with angular 
frequency ujp and power Pp, is injected into the feedline. The 
reflected signal is split and measured by a spectrum analyzer 
in the frequency domain, and by a lock-in amplifier (LIA). 
The LIA is tuned to the frequency of the slow magnetic field 
modulation 6^ cos(ajLiAt) that is applied through a wire. The 
inset shows an electron micrograph of the microbridge. 



a function of the input frequency. Near a magnetic field of 
Bf^^ = 0.215 T the resonance frequency of the resonator 
was shifted and the Q- factor decreased (see Fig. The 
results obtained by us are similar to those presented by 
others dealing in the coupling between a paraiiiagnetic 
sample and a linear superconducting resonator |14Hl6l |. 

We then examined the response of the resonator at 
the non linear regime. We excited the resonator by in- 
jecting a monochromatic tone at frequency Wp which is 
close to the resonance frequency, and power in the vicin- 
ity of the modulation threshold. The static magnetic field 
was slowly ramped around Bf^^, while a low frequency 
(1.23 kHz) low amplitude (20pA) ac-current was injected 
into the wire (as in conventional CW induction detec- 
tion schemes) providing the magnetic field modulation. 
A lock-in amplifier (LIA) was tuned to the modulation 
frequency and measured the envelope of the signal re- 
flected off the resonator. Slightly below the modulation 
threshold, the response was linear and power reflected 
off the resonator was proportional to the 6*11 parame- 
ter (the reflection coefficient); therefore the LIA mea- 
sured the derivative of the Sn parameter with respect to 
the magnetic field. The ESR-induced change in the LIA 
signal was seen at the same value of Bf^^ as in direct 
Sii measurement (Fig. ^h, blue line). When the input 
power was set exactly at the modulation threshold, the 
LIA signal increased significantly by a factor of up to 
100 (Fig[4]-b, green line). As we scanned the static mag- 
netic field we found that the non linearity power thresh- 
old shifted; at Sf^R ^j^g 

power threshold increased by 
~ 1.2 dB relative to the case where the static field is far 
from the resonance value. Thus, by scanning over the Pp 
- Bz plane, and measuring the LIA signal we can obtain 
the ESR spectrum, (see Fig. |3]-a). 

To analyze the results, we consider an ESR sample 
placed near a current anti node in a stripline MW res- 
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Figure 2. (color online) Panel A: Self modulation stability di- 
agram. The stability map, which is obtained by Eq. Q con- 
tains 4 stability zones: the superconducting mono stable zone, 
MS(SC), the normal conducting mono stable zone, MS(NC), 
the bi stable zone, BS, and the astable zone, AS. We com- 
pare the case were B- is far from Bf^^ (blue-solid lines) and 
the case of B^ — Bf^^ (red-dashed lines). The experiment 
in Fig. 13] is done along the black line, at input power range 
— 3dBm < Pp < — IdBm, and slightly above the resonance 
frequency. Panels b-c: Data from spectrum analyzer at the 
monostable (C, Pp = — 3dBm) and astable (b, Pp = —1 dBm) 
zones, when Bz is far from B?^^. 
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onator embedded with a microbridge. An external static 
magnetic field tunes the Larmor frequency of the ESR 
sample close to the resonator's resonance frequency. The 
dynamics of the resonator in such a case can be captured 
by two coupled equations of motion. The full derivation 
is provided in [l^. Below we describe them in brief. 

The resonator is driven by a coherent tone a™e-*'^p* 
injected into the test port, which is weakly coupled to 
the resonator, where a*" is a constant amplitude and Wp 
is the drive angular frequency. The mode amplitude A 
can be written as A = Ae"'"^?*, where A (t) is a complex 
amplitude, which is assumed to vary slowly on the time 
scale of l/ojp. In this approximation, the equation of 
motion of A reads 



dA 
dt 



[i (wp - cjo) - 7] ^ - *i 



where luq is the resonance frequency, 7 ~ 71+72+73, 7i is 
the coupling constant between the resonator and the feed- 
line, and 72 + 73 is the damping rate of the mode, where 
72 denotes the dumping rate of the microbridge, and 73 is 
the dumping rate of all other loss mechanisms. The res- 
onator's Q-factor is defined as Q — uJof-f. The term c'" 
represents a random-phase input noise with correlation 
function: (c*") = and (c*"(Oc™*(t')) = 27^(5 {t - t'), 
for the case of thermal equilibrium at high temperature 
[kBT ^ fiwo), where /cb is Boltzmann's constant. 

We consider the case where non linearity is generated 
by a local hot-spot in the resonator, i.e the microbridge. 
The microbridge is assumed to be sufficiently small, so 
that its temperature T may be considered homogeneous. 
The temperature of all other parts of the resonator is 
assumed be equal to that of the coolant Tg. The heat 
balance equation for the microbridge reads 



C— ^-H{T- To) + 2r7^p72 



A 



(2) 



where C is the thermal heat capacity and H is the heat 
transfer coefficient. 

The coupling mechanism between Eq. ((Ij and Eq. ^ 
is based on the resonator parameters' dependence on the 
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Figure 3. (color online) The resonator's resonance frequency 
and Q factor extracted from Sn measurements under scan. 
Experimental results (solid blue line) are compared with the- 
oretical prediction based on Eq. ^ (dashed red line). The 
parameters used in Eg. @ are: ti = t2 — l/61MHz/7r, 
iV = 2 • lO^^spins/m^ Tesr = 4.2 K, and = 8.6 ■ 10"^ 
The simulation takes into consideration a single homoge- 
nous broadened ESR line, while in experimental data g-factor 
anisotropy and hyperfine effects are observed [l^ . 



impedance of the microbridge (T9| . which is dependent 
on its phase, i.e., superconductive (SC) vs. normal con- 
ductive (NC). We assume the simplest case, where the 
dependence is a step-function at the critical temperature 
T = Tc, namely wo,72 take values a;o,sj72,s for T < Tc 
and values wo,n, 72, n 

for T > Tc [Hi. ' 



In general, the two coupled Eqs. ([T]) and ©, have 
two attractor sets [l^: {^cxj.s, Too.s} (the "super attrac- 
tor") and {Ago, n, Too, n} (the "normal attractor"). When 
T < Tc the system evolves toward the super attractor, 
whereas, when T > Tc it evolves toward the normal at- 
tractor. The stability of the super (normal) attractor is 
dependent on the relation Ttx),s < Tc (Too.n > Tc)- Conse- 
quently, four stability zones can be identified in the plane 

I vs. pump frequency Wp 
(see Fig. [5]). In the monostable zones, either the SC 
phase or the NC phase is locally stable; in the bistable 
zones, both phases are locally stable. In the astable zone, 
on the other hand, none of the phases are locally stable, 
and the resonator oscillates between these two phases. As 
the two phases significantly differ in their reflection coef- 
ficients, the oscillations are translated into an amplitude 
modulation of the reflected pump tone. The borderlines 
between the four stability zones are given by: 



+ c*", (1) of pump power Pp = huj-p 
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Figure 4. (color online) ESR spectrum obtained in the linear 
and non-linear regimes. (A) The signal reflected off the res- 
onator, as a function of the static magnetic field and the input 
power, measured using a LIA at 1.23 kHz. The ESR spectrum 
can be obtained from the non-linear response readings (local 
maxima of the graph), and is compared to analytical results 
obtained from Eqs. pa|) and (O (dashed white line). The 
fitting parameters are H = and Tc — To = 70 K. 

(B) Cuts from the two dotted lines which are marked in the 
colormap: — 2.8dBm (solid-blue) and — 2.2dBm (solid-green) 
input power, showing linear and non linear response respec- 
tively. Note the discontinuity in the vertical axis, which was 
created in order to make visible the relatively small change in 
the linear regime. 
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H{T, - To) 
47i72,s 



^0,sf + Is] , 



(3a) 



47172,11 

Eq. pap determines the border line between the super 
conducting mono stable zone and the astable zone which 
is the modulation threshold. The experiments were done 
in the vicinity of this borderline (Fig. ^a). 

We now consider an ESR sample under static external 
magnetic field = B^z and MW magnetic field B^; = 
Bx cos{ujpt)x, such that B^ ^ B^- Under such conditions 
the susceptibility of the ESR sample is given by 0, [l3| : 



X = 

X = 



(Wp - uj{)tl 



1 

'2Ynujp-uj,)^q + lls^Blht2 

1 t2 

2 1 + (wp - wO'^i + 7isR^'^ii2 



t^iXo, (4a) 
^\Xo, (4b) 



where uj\ = 7esr52 is the Larmor frequency, ti and t2 
are the ESR longitudinal and horizontal relaxation times 
respectively, and xo is the static susceptibility which is 
given by xo = Mo-^7esr^^/4^b'7esr for the case of spin 
1/2 system with g-factor 2, where Tesr is the ESR sample 
temperature. The interaction with the spins gives rise 
to changes in the resonator's resonance frequency, wq, 
the dumping coefficient, 7, and the quality factor, Q, 
according to the expressions [l|: 



AQ/Qo « -A7/7 « - 



ioVX 



(5a) 
(5b) 



where ry is the filling factor of the sample In the lin- 
ear regime when analytical results are compared with the 
experimental results in Fig. [3] good agreement is found. 
As can be seen from Eq. ([5a|) the ESR shifts the border- 
line between the mono-stable and astable zones. Again, 
when comparing the analytical results with the experi- 
mental results in Fig. |4l good agreement is obtained. 

A potential advantage of our methodology can be 
explained by comparing it to the conventional linear 



method. In ESR inductive detection with a resonator 
having a linear response the output signal is relatively 
small and consequently various readout elements such as 
cryogenic amplifiers, LIAs, etc. are commonly employed. 
Reducing the noise of such readout elements to cryogenic 
temperatures is a major engineering challenge [2l|, [2^ . 
In our scheme, on the other hand, the ESR signal am- 
plification is achieved by the resonator's intrinsic behav- 
ior, and consequently no further active amplification is 
needed. A second potential advantage of our methodol- 
ogy is important for the case of paramagnetic materials 
with broad resonance lines. This advantage can be un- 
derstood in relation to LIA measurement and field modu- 
lation. Commonly, LIA measurement is applied in order 
to reduce 1// noise. In the present experiment the LIA 
detection is carried out with respect to a modulation in 
the static magnetic field with low frequency. This mod- 
ulation brings the resonator in/out of the nonlinearity 
boundary. Such detection scheme is the most common 
in CW ESR and facilitates its high sensitivity. How- 
ever, for paramagnetic materials having broad resonance 
lines reaching 1000 G P, 01 and more, the conventional 
field modulation at a large amplitude is not possible, as 
it causes excessive heat and mechanical vibrations. This 
problem can be solved in our scheme by employing ampli- 
tude modulation of the driving MW frequency (instead of 
field modulation), for which the width of the ESR spec- 
trum does not limit the detection (see Fig. [5]). 

In summary, we have introduced a new methodology 
for induction detection of ESR utilizing the non linear 
response of the MW resonator. We have presented the 
phenomenon experimentally and theoretically. Our ini- 
tial setup is not optimal for ESR measurements of small 
samples, but in the future, we plan on combining this 
method with compact MW resonators whose magnetic 
field is confined to a small volume 
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